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' We study the many-body effects in a two-subband quasi-two-dimensional electron system in a 

quantizing magnetic field at filling factor three. A manifestation of these effects in the phonon 
^ ' absorption spectroscopy is discussed. The electron system is mapped onto a two-level system with 

Q , the separation between levels determined by the intersubband splitting and the cyclotron energy. 

■ The electron-electron interaction enhances the excitation gap, which exists at all values of the 

interlevel splitting. This results in a single-peak structure of the phonon absorption rate as a 
, function of magnetic field, instead of the double-peak structure for non-interacting electrons. 
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I. INTRODUCTION 



The behavior of a two-dimensional (2D) electron system in a quantizing magnetic field is completely determined by 
the interaction part of the Hamiltonian. This results in several interesting phenomena in such a system, for example, 
, including the Fractional and Integer Quantum Hall Effects. □ A wide range of different experimental techniques kas 
been used to study this system and its excitations. One of these methods is the acoustic-phonon spectroscopy.H^El 
I \ The specific feature of the electron-phonon interaction in the 2D system in a magnetic field is a strong suppression of 
■ O • phonon induced transitions with the phonon momentum larger than the inverse magnetic length, go ~ 1/1^ Due to 
^ I the energy conservation the singic-phonon transitions are only possible if E{qo) « hsqo, where s is the speed of sound 
and E{q) is the dispersion of the neutral excitations of the electron system. By changing E{qo) one can switch off or 
on acoustic phonon effects in the magnetically-quantized 2D system. 

At certain conditions a sharp magnetic field dependence of the electron-phonon interaction can be achieved in a 
two-subband system. When the energy splitting between the first and the second subbands becomes close to the 
cyclotron energy the second Landau level of the first subband has almost the same energy as the first Landau level of 
the second subband. An interesting feature of this system is that the separation between these close Landau levels, 
A, can be changed by changing the magnetic field. Under these conditions if there are electrons in the second Landau 
level of the first subband, i.e. the electron filling factor is greater then two, the phonon-induced effects become strong 
when A « hs/l. This results in a double-peak structure of the dissipative electron conductivity as a function of 
A or as a function of magnetic field.El The increase of the electron density in this system increases the role of the 
electron-electron interaction which results in the renormalization of the excitation gap A and has a considerable effect 
on the electron-phonon coupling. In what follows we analyze the case when the total electron filling factor is equal to 
three. We show that in this case the excitation spectrum has a gap at all values of interlevel splitting A. As a result 
C [ the acoustic-phonon absorption of the two-subband system has a single peak as a function of the magnetic field. By 
• lowering the electron density one should observe the transformation of a single-peak dependence into a double-peak 
' one, which corresponds to a non-interacting case. 
O ■ We consider below the quasi-2D electron system confined in a single heterojunction. The effects related to the 
J-^ ' intersection of different Landau levels belonging to different subbands were studied both theoretically and experimen- 
tally for double quantum well structures.cl Unlike the double-well case the intersection of levels in the two-subband 
single heterojunction occurs at much higher magnetic field, which makes all field-dependent many-body effects more 
pronounced. 



II. GROUND STATE AND EXCITATIONS 



At filling factor vq = 3 the electrons completely occupy the lowest Landau level (with both spin directions) of 
the first subband. The only effect of the electrons in this level on the electrons in the higher levels is the exchange 
renormalization of the interlevel splitting A. Because this renormalization does not depend on the momentum of the 
electron in the next levels in what follows we do not take the background (completely occupied) level into account and 
consider A as the value, which already includes the renormalization due to the background electrons. As a result the 
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system becomes the two-level electron system with the effective filling factor v = 1 and the splitting between levels 
A, see Fig.l. These levels are the second Landau level (n = 1) of the first subband and the first Landau level (n = 0) 
of the second subband. Both levels have the same spin, e.g., Sz = +1/2 for a negative ^-factor. 

It should be mentioned that the mapping of the i^q — 3 system onto the two-level problem is valid for any value of 
Zeeman splitting. Due to exchange interaction the ground state is a mixture of the states with the same spin only. 
Since we are interested in the phonon induced transitions, which conserve the electron spin, we should consider the 
excitations without spin reversal. As a result only the electron states with the same spin are relevant to the considered 
problem. Because both levels have the same spin we can treat electrons as spinless particles. For convenience the levels 
are numbered by the corresponding subband numbers, = 1 or 2. To simplify the expressions the energy is measured 
in the Coulomb units ec = / kI and the length is in units of magnetic length /. Following the standard procedureB 
of rewriting the interaction part of the Hamiltonian in the momentum representation, the total Hamiltonian of the 
two-level electron system can be written in the form: 

+ 2 ^ ^ ('^iA'4M/^2) 2^ '-'fei+9„,A.i'-'fe2,M2'-''=2+9«^A'3<-fel,A'4 : UJ 

{m} 9^.9!/ fci.fc2 

where the Landau level number is equal to rij = 1 for [li — \ and = for /i^ = 2, i.e. = 2 — p^i. We use the 
Landau gauge with vector potential A — (0, Bx, 0) and Cj^^^, C^^^ are the creation and annihilation operators of the 
electron in the state ipn,ky,fj.- 

Ipn.kyA^^ViZ) = Xf^{z)—^(t>n{x - ky) , (2) 

where Xfj.{z) is the envelope function of the /ith subband; (j)n{x) is the nth harmonic oscillator function. 
In Eq. (1) we introduced the notations: 

^i7™f('?) = J^mM4P3M2('z)G„,„,((?)G„3„2(r) exp(^-0 , (3) 

wherei 

with q = qx + iqy, q — \q\ and being a generalized Laguerre polynomial. The only functions relevant to our 
problem are 

f <f'\ q ( q^ 

Go,o(g) = exp 1^-— j , Go,i((7) = ;^exp ^^- — 
Gl,l('Z)= (l-y)exp^ 

The modification of the Coulomb interaction due to the finite extent of the electron wave functions in z-direction is 
given by the expression: 

/•oo /"OO 

■F)«lA'4M2M3(9) = / / rfziC?2:2e~'^l'^^~^'lxMl(^l)XM4(^l)XM2(2:2)XM3(^2) ■ 



We use the Fang-Howard approximation0 with the parameter h for the envelope functions of the first and second 
subbands: 

X2W = /f^(^-^)exp(-i6z 
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where the function X^iz) is found from the condition that Xi(z) and X'liz) are orthogonal. 

The first problem we are interested in is the ground state of the system. It is clear that if the splitting A between 
the layers is larger then the characteristic inter-electron Coulomb interaction (A > 1) then the ground state is a 
completely occupied level with /i = 1. Similarly, if A < — 1, the ground state is a completely occupied level with 
/i = 2. With changing A the transition between these two phases occurs. This transition can be either sharp or 
smooth and in the last case the transition occur through the new intermediate phase. The example of sharp phase 
transition in the two-level electron system with the total filling factor = 1 is given in Ref. ^ where the two levels 
correspond to different Landau levels of the same subband but with opposite spin. There is aiso an example of a 
smooth phase transition in the two- level system of composite fermions with filling factor i/ = l.tj The two levels in 
Ref. |l^ are the same as the levels in Ref. The only difference from Ref. |l^is in the type of interaction between 
the particles. Our system is different from the systems considered in Refs. in two aspects: the levels in our 

case correspond to different subbands and the spin for both levels is the same, i.e., we have an additional exchange 
interaction between the electrons in the different levels. 

To find the ground state of the system in the intermediate range of A we follow the method of Ref. n2 and use 



the standard Hartree-Fock approximation. We assume the non-zero average values of 



^k.2^k,2) and 



i^k+Q,2^ which do not depend on k. In general, the parameter Q which minimizes the Hartree-Fock Hamiltonian 

is not equal to zero. For each value of k we introduce the new wavefunctions which are the eigenfunctions of the 
Hartrcc-Fock Hamiltonian. The creation and annihilation operators corresponding to these functions are a^^, ak,i, 
j = 1, 2. The new functions are related to the original ones by the rotation by angle 9: 



Ck,i = afc,i cos(6'/2) -I- afc,2 sin(6'/2) , 
Ck+Q,2 = -flfea sin(6'/2) + 0^,2 cos(6'/2) 



(5) 



Assuming that ak^i corresponds to the lowest energy state we write the average: \CLj^iO,k,jj — SuSji. Substituting 
expression (5) into the Hamiltonian (1) we get the Hartree-Fock Hamiltonian as a function of 9 and Q: 

1 



H 



HF 



[-A - €\\ cos^ {912) - sin^ {9/2)] cos^ {9/2) 



^22 



sin^ (61/2) 



-12 



cos' 



{9/2)] sin^ {9/2) 



-sin2 (0/2)cos2 {9/2)Vi2{Q) 



where ej^*^^ is the exchange energy of an electron in the iith Landau level of the jith subband interacting with the 
electrons of the same spin in the filled i2th Landau level of the j2th subband: 
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dkFj^j^j-^^j^{k) Gi-^i2{k) 



exp 1 - 



(6) 



The potential Vi2{Q) is given by the expression: 



Vi2{Q) = / dk 



Fii22{k)Gn{k)G„o{k)Jo{kQ) - Fi2i2{k) Gi2{k) J2{kQ) 



expl-y 



Analyzing the Hartree-Fock Hamiltonian, one can show that the energy minimum can be achieved only at = 
or = TT, which corresponds to a completely occupied level /i 1 or /i = 2, respectively. This means that there are 
only two different phases with the a sharp transition between them. Comparing the Hartree-Fock energy at = 
and 9 — n, we obtain the inter-level splitting A* at which the transition between two phases occurs: 



(7) 



For A > A* the ground state is the completely occupied level with /i = 1 and for A < A* the ground state is the 
completely occupied level with [i = 2. 

To find the energy spectrum of neutral excitations we follow the method proposed in Ref. ^ for a double layer 
system at electron filling factor equal to 1. We introduce the normalized cxciton operator, which creates an electron 
in the empty level and a hole in the filled level: 



-4+ (9) 



(8) 
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where q = {qx, Qy) is the momentum of the electron-hole pair, A'' is the number of electrons and we assume that the 
level with = 1 is filled. Then we calculate the energy spectrum of such excitations as an average of the Hamiltonian 
(1) in the state |0); 

E{q) = {A{q)HA+{q)) = A + e}} - el\ + |Fi2i2(g)cxp {-q' 12) 
- dk{l-^^ Jo(fcg)^ii22(fc)exp(-fcV2) 
= (A - A*) + i (6^° + e\\) - ell + f ^1212(9) exp {-q^ /2) 

-J^ Jo(fc(7)^ii22(fc)exp(-fcV2) , (9) 



where we rewrite the E{q) in terms of the difference (A — A*). Because the Hamiltonian (1) does not conserve the 
number of excitons we should take into account the processes which involve several excitons. Such processes are 
important at small pSialues of the level splitting A. To take into account the multi-exciton states we approximate the 
Hamiltonian (1) bjO 

H^Y. [Eiq)A+iq)A{q) + i [X{q)A{q)A{~q) + X* iq)A+ {q)A+ {-q)] | . (10) 

The function X{q) is defined by the equation: 

X{q) = {HA+{q)A+{-q)) = e^^^ ^F.^M exp (-gV2) 

-e^^^y^ (ifc|-J2(Wi2i2(fc)exp(-fcV2) , (11) 

where <f> — aTctan{qy / q^) ■ After diagonalizing the Hamiltonian (10) using the Bogoliubov transformationO we obtain 
the energy spectrum of the elementary neutral excitations: 



Eoiq) = VE'i'l)-\m\' ■ (12) 

If the level /i = 2 is occupied in the ground state, i.e. A < A*, then the Eqs. (11) and (12) are the same and the 
only difference is in Eq. (8) where we should write —(A — A*) = |A — A*| instead of (A — A*). 

The numerical calculations show that the lowest energy excitation (Eq. (12)) has a zero momentum. The energy 
of this excitation can be considered as a renormalized splitting, A^, between the levels of the double- level electron 
system. It follows from Eqs. (9), (11) and (12) that 

A, = ii;o(0) = |A - A*| + i (e™ -I- ell) - 41 - dk [l ~ ^) F.^k) exp {-^ /2) . (13) 

The splitting A^ reaches its lowest value at A = A*. The renormalization of the excitation energy in a 2D elec- 
tron system (^pfectpr enhancement) due to the inter-electron interaction was well studied both theoretically and 
experimentally.E3~E3 At the same time there are excitations of electron systems which are not altered by the electron- 
electron interaction. For example, Kohn theoremt^ tells us that the energy of a cyclotron excitation does not depend 
on the electron density. 

In Fig. 2(a) the critical level splitting A* is shown as a function of the Fang-Howard parameter 6 by a solid line. 
With increasing b the spreading of the electron wave function in z direction becomes smaller. This results in increasing 
A*, which can be understood from the Eqs. (6) and (7). The first term in equation (6), e™, is determined by inter- 
electron interaction in the first Landau level of the second subband, while the second term e}} is determined by the 
inter-electron interaction in the second Landau level of the first subband. Since the spreading of the electrons in 
the second Landau level is broader than the spreading of the electrons in the first Landau level, the 2D form-factor 
of the inter-electron interaction, which is given by Gi^^i^{q) (Eq. (4)), is larger for the electrons in the first Landau 
level. Similarly, due to a different spreading of the wave functions in z-direction for the electrons in the different 
subbands the form-factor F^-^^^^^^^^{q) is larger for the electron in the first subband. The competition between these 
two tendencies gives rise to the increase of A* with increasing b. We can see from Fig. 2a that at b k, 1.6 the effects 
of the different spreading in ^-direction and in the (a;, ?/)-plane cancel each other, so that A* — 0. 
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The dotted hne in Fig. 2 (a) shows the dependence of the renormahzed sphtting on the parameter b for A = A*. 
We can see again that with decreasing of the spreading of electron wave function in z direction the splitting A^ 
increases. This increase of A^ has the same reason as the increase of A*. 

In Fig. 2(b) the dispersion curve -Eo(<z) (Eq. (12)) is shown for A = A* and for different values of parameter b: 
6=1 (solid line) and b — 2 (dotted line). We can see again that the excitation energy increases with increasing b. 
The effect becomes stronger for larger values of momentum q. 

In Fig. 2(c) the dispersion curve Eo{q) is shown for 6 = 1 and for two different values of A: for A = A* (solid line) 
and for A = A* + 0.05 (dotted line). The increase of A has no significant effect on the form of the dispersion curve 
due to a small value of the function X{q}. As a result the dispersion curve Eolk) has an almost linear dependence on 
splitting A for all values of momentum q and can be approximated by Eo{q) « E{q). 



III. PHONON SPECTROSCOPY 



In this section we study the interaction of the quasi-2D double subband electron system, discussed in the previous 
section, with acoustic phonons. We consider the case of the absorptionjDhonon spectroscopy when the 2D electrons 
are subjected to an external beam of non-equilibrium acoustic phonons Eij 

Below we are using the isotropic Debye approximation in which the phonon frequency has a linear dependence on 
the wave vector: 

ujj{K) ~ SjK , 

where Sj is the speed of sound, j is labeling the phonon modes, j = 1 for longitudinal mode and j = 2, 3 for two 
transverse modes. The wave vector K is the three dimensional (3D) vector, K = \K\. For convenience we label the 
3D vectors by the capital letters and their projections by the corresponding small letters, K = {k,kz)- 

The electron-phonon interaction Hamiltonian can be written in terms of cxciton operators A^,A (Eq. (8)): 



e—ph 



E 



M,iQ) 



Z{qz) A{q)d+{Q)+A+{q)d,{Q) 



(14) 



where is the creation operator of a phonon in the jth mode, y is a normalization volume; Mj(Q) are the matrix 
elements of electron-phonon interaction which are determined by the deformation potential and piezoelectric coupling. 
In GaAs the matrix elements Mj{Q) are usually written in the formt2l: 



2/90 



-f3 



(15) 



where po js the mass density of GaAs, f3 and are the parameters of piezoelectric and deformation potential 
couplingsj3 respectively; is the polarization vector of the jth phonon mode. 

The form-factor Z{qz) in Eq. (14) is determined by the electron spreading in z direction: 



Ziqz) = / dze'''^'xiiz)x2iz) 



(16) 



To illustrate the specific feature of the electron-phonon interaction in our two-subband system we study the total 
rate of absorption of phonons from an external pulse by the 2D electrons. At low enough temperature the rate of 
absorption is given byll3 



27r 



E 



dO ^ ^ 2 

-^S{Eo{q)~sQ)n,{Q) M,{Q)Z{qz) Roi{q) , 



(27r 



(17) 



where nj{Q) is the phonon distribution function. The function i?oi('?) corresponding to the transition from the first 
to the second Landau level is given by 



i?oi(<z) = IGoAmW 



^,-,V2 

2 



(18) 
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In Eq. (17) we assume that the energy relaxation time of the 2D electron system is much shorter than the charac- 
teristic time of the phonon absorption, which as shown below is of the order of 10^^° s. 

From the expression for Roi{q) (Eq. (18)) one can see that the phonon absorption occurs mostly for the 2D 
momentum q being in the range of g ~ 1 2. In this range of q the energy spectrum Eo^q) is almost dispcrsionless 
(Fig. 2(b)), i?o(<?) ~ 0.08 + |A — A*| for 6=1. Energy conservation requires the energy of absorbed phonon to be 
equal to E^lq), i.e. S\/ q^ + q^ — Eolq). Taking into account that the phonon momentum q^ is restricted by the 
parameter b {\qz\ < b) and that the speed of sound is about 0.03 for longitudinal and 0.02 for transverse phonons, we 
see that the phonon absorption is suppressed in the two-subband system. This suppression is completely due to the 
many-body effects, since the many-body gap, EQ{q)ec, is greater than hs/l. 

Let us assume that the non-equilibrium phonon distribution function is isotropic. Then, substituting Eqs. (15) and 
(18) into the Eq. (17) and performing the integration, we get: 

1 _ Wafcs 

Tabs n{Eo{qo)) 

9 M, 



- -J exp(-9?(l-«2)/2) 

8 (Zi 



% / du{l-u^)\Z{uq,f 
^1 Jo 



L 

(?3 Mp 



+ / du{l~u'y{9u^ -2u' + l)\Ziuq2)\exp{-qi{l-u')/2) , (19) 

S2 8 Jo 

where qi — Eo^qo)/ si, q2 = £'o('?o)/'S2 and go is taken on the plateau of the dispersion curve Eo{q) (Fig. 2(b)). Here 
the constants Md and Mp, characterizing the deformation and piezoelectric couplings, are: 

Md = 
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In Fig. 3 the absorption rate is shown as a function of A — A* for the two values of parameter b: b = 1 (solid line) and 
b = 1.5 (dotted line). Both curves have a single-peak structure and show a weak dependence on b. The energy Eo{qo) 
is increasing with increasing b (Fig. 2(b)). This results in the decrease of the phonon absorption rate. However, the 
increase of b is favorable for the absorption of phonons with larger g^, because the form- factor Z{qz) effectively cuts 
off the phonon momentum g^ at g^ = 6. These two factors counteract each other, which results in a weak absorption 
rate dependence on parameter b. For the focused phonon beam, the first effect becomes more important and the 
decrease of the absorption rate with increasing b should be observed. 

The rate of the phonon absorption for the non-interaction electron system is shown in Fig. 3 by the dashed line. 
The main difference between the interacting and non-interacting systems is in the shape of the dependencies of the 
absorption rate on the level splitting. These dependencies have the double-peak and single-peak structures for non- 
interacting and interacting systems, respectively. For the non-interaction system the excitation spectrum is just the 
bare gap |A| and the excitation energy has any value starting from 0. The spectrum becomes gapless at A = 0. In 
this case the maximum of the absorptioa,rate corresponds to the condition: |A| ~ s/l. This results in a double-peak 
structure for the non-interaction system.Q For the interacting system there is an energy gap for any value of A. Thus, 
the absorption rate as a function of A has a single-peak structure. 

If the electron filling factor v is changed from = to = 1 one should expect the transformation of the double-peak 
structure into a single-peak one, because at small filling factor the inter-electron interaction becomes less important. 
This transformation occurs as the merging of the two peaks into a single peak. Assuming that the interaction gap 
of the 2D electron system is proportional to i^, one can see that with increasing v the separation between peaks 
decreases linearly, and the two peaks transform into a single peak a.t f Ki 0.75 for 5=1. The strength of the peaks 
increases linearly with increasing v, because the phonon absorption rate is proportional to the electron density. The 
transformation of the double-peak into a single-peak structure with increasing electron filling factor can be observed 
experimentally. The splitting between the levels A is the difference between the intersubband splitting and cyclotron 
energy with some corrections due to exchange interaction with the background electrons. Therefore, this splitting can 
be changed by changing magnetic field. 

At the same time the gate voltage should be adjusted to keep the electron filling factor constant. However, even for 
the changing filling factor at fixed electron density one should expect the single-peak structure of the phonon absorption 
curve as a function of magnetic field if the electron filling factor is close to ly = 1. The double-peak structure appears 
when the electron filling factor is smaller then ^ 0.7. For a constant electron density the single-peak and double-peak 
structures are both asymmetric when the difference A — A* changes sign. 
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IV. SUMMARY 



We have studied a two-subband electron system in a quantizing magnetic field at the total filling factor uq = 3. 
When the second Landau level of the first subband is close in energy to the first Landau level of the second subband 
the system becomes effectively a two-level electron system with filling factor v = 1. Depending on the value of splitting 
between the two levels the system can be in one of the two possible ground states: for A > A* all the electrons are 
in the first subband and for A < A* all the electrons are in the second subband. The inter-electron interaction 
renormalizes the critical splitting A* which should be zero for a non-interacting system. Another effect related to 
the electron-electron interaction is a non-zero gap for any value of level splitting A. This effect should be important 
for the phonon spectroscopy experiments. We show that due to the non-zero many-body gap the phonon absorption 
rate has a single-peak structure as a function of the splitting between the levels or as a function of magnetic field. 
At small electron filling factor i^, when the system can be considered as a non-interacting system, the double-peak 
structure is expected. The transformation of the double-peak structure into a single-peak one occurs when the filling 
factor is changed from u = to u = 1. Our calculations show that the two peaks merge at filling factor u 0.7. 
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FIG. 1. Single-electron energy levels of the two-subband system. Each level is numerated by the pair of numbers: the Landau 
level number and the subband number. The arrows indicate spin directions. The interlevel splitting is A. Zeeman splitting is 
A.. 
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FIG. 2. (a) The critical level splitting A* (Eq. (7)) and the renormalized splitting A,, at A = A* (Eq. (13)) are shown 

as the functions of the Fang-Howard parameter b by solid and dotted lines, respectively, (b) The dispersion curve Eo{q) (Eq. 
(12)) is shown for A = A*. The numbers near the lines are the values of the parameter b. (c) The dispersion curve Eo{q) (Eq. 
(11)) is shown for 6 = 1 and different A: A = A* (solid line) and A = A* -|- 0.05 (dotted line). Both q and b are in units of 
inverse magnetic length, l/l; the energy is in units of ec- 
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FIG. 3. The phonon absorption rate wo = l/^abs (Eq. (19)) as a function of A — A* for 6 = 1 (solid line) and b = 1.5 (dotted 
line). The absorption rate for the non-interacting system is shown by the dashed line. The absorption rate u!o is in units of 
10^° s"\ A - A* is in units of ec- 
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